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ABSTRACT 



This paper presents a multicolor study of the nearby cluster of galaxies Abell 168 
[z = 0.045) with 13 intermediate-band filters in the Beijing- Arizona- Taiwan- 
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0> ! Connecticut (BATC) filter system. After a cross-identification between the pho- 

\Q '. tometric data obtained from the BATC and the Sloan Digital Sky Survey (SDSS), 

a^ 

O . 

. eludes 121 spectroscopically confirmed member galaxies. The technique of photo- 

. metric redshift has been applied to all these galaxies with combined 18-band (13 

p i' from BATC and 5 from SDSS) spectral energy distributions (SEDs), in order to 

Q ■ perform a faint membership selection in Abell 168. As a result, 255 galaxies are 

• newly selected as the member candidates. 

; Based on the enlarged sample of cluster galaxies, the spatial distribution and 

dynamics of Abell 168 are investigated. In the light of the spatial distribution of 
^ ' the member galaxies and the 0.2 — 3.5 keV X-ray image by Einstein observatory, 

I it seems that Abell 168 consists of two merging subclusters with a relative radial 

velocity of 264 ± 142 km s~^. With the help of ROSTAT software, a detailed inves- 
tigation of the dynamics shows the intrinsic difference in the velocity distributions 
for these two subclusters. Under a linear two-body model, they are found to be 
a gravitationally bound system with 92% probability. The slight deviation of the 
local velocity distribution from the overall distribution in the central region may 
suggest a picture of head-on collision in the projection plane. By using a double 
Schechter function, the close examination of the luminosity functions reveals the 
different galaxy content of these two subclusters, which implies the merging process 
is possibly at the early stage. 



^National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, P. R. China 
^Department of Physics, Nanjing Normal University, NingHai Road 122, Nanjing, 210097, P. R. China 



- 2 - 



Subject headings: galaxies: clusters: individual (Abell 168) — galaxies: distances 
and redshifts — galaxies: kinematics and dynamics — galaxies: luminosity function, 
mass function 

1. INTRODUCTION 

The increasing volume of observational data for galaxy clusters (GCs) in multiple wave- 
bands provides important constraints upon the formation of large scale structures and cosmol- 
ogy (Bahcall 1988; Dressier & Gunn 1988; Schindler 2001). The investigations on the luminosity 
function and the spectral features of member galaxies in a rich cluster will help us to understand 
not only the structures and dynamics of the whole galaxy cluster but also the formation and 
evolution of galaxies under such a dense environment (Dressier 1984; Poggianti 2002). 

Spectroscopic observation in the optical band is still the most straightforward and powerful 
approach to the membership determination. However, the reliable spectroscopics for the faint 
galaxies in a cluster is still an unattainable goal for the large telescopes. With the development 
of photometric redshift technique, the redshifts of faint galaxies can be significantly determined 
according to their SEDs. The BATC sky survey is specially designed for this purpose. In 
this paper we present a large-field multicolor study of the galaxy cluster Abell 168 based on 
the BATC observations with 13 intermediate-band filters, which covers the whole optical band. 
Improvement of the accuracy in photometric redshift estimate with our multicolor system allows 
us to study the spatial and the dynamic properties of the nearby clusters of galaxies (Yuan et 
al. 2001; Yuan et al. 2003). 

Abell 168 is a nearby {z = 0.045) cluster of galaxies located at 01^15'^10^2 +00°12'18" 
(J2000.0) with richness of I— HI (BM classification). Its distance module is 36.45 by assuming 
Qo — 0.5 and Hq — 75kms~^Mpc~^ (also used throughout this paper). The luminosity func- 
tion of Abell 168 was remarked as the best example of non-Schechter form by Dressier (1978), 
which was confirmed by Oegerle, Ernst & Hoessel (1986). Based on the positional discrepancy 
of ~ 10' (i.e., ~ 540 kpc) between the X-ray and the optical center, Ulmer, Wirth & Kowalski 
(1992) suggested that Abell 168 was formed by the collision of two clusters with approximately 
equivalent sizes. It is reported that the X-ray emission of Abell 168 has two peaks with sur- 
face brightness of 9x10^^"^ and 3x 10^^"^ erg s~^cm^^, respectively (Tomita et al. 1996). They 
supposed that the star-forming processes will be triggered when gas-rich galaxies rush into the 
higher density regions. However, Tomita et al. (1996) failed to find the enhanced fraction of 
blue galaxies in the region between these two X-ray peaks. 

This paper aims at not only selecting faint member galaxies of Abell 168 by photometric 
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redshift technique, but also investigating the substructures in the spatial distribution and the 
line-of-sight velocity distributions. We wish to provide the direct evidence of the ongoing merge 
event of Abell 168 from our analysis of multicolor photometries. This paper is structured as 
follows: In §2, we describe the BATC observations, the data reduction. The technique of 
photometric redshift is used to estimate the redshift for all galaxies, then the membership 
selection of Abell 168 is performed. In §3, the substructures of Abell 168 are explored on the 
basis of the member galaxy distribution and the X-ray emission of hot gas. The dynamics of 
substructures is also studied. In §4 the luminosity functions (LFs) of Abell 168 are analyzed. 
The summary is given in §5. 



2. DATA 

2.1. BATC Observation and Data Reduction 

The BATC program is a large- field (58' x 58') multicolor sky survey based on the 60/90 cm 
f/3 Schmidt Telescope of National Astronomical Observatories, Chinese Academy of Sciences 
(NAOC). A Ford Aerospace 2048x2048 CCD camera with the pixel size of 15 //m ( T'.T/pixel) 
is mounted on the telescope. The BATC filter system contains 15 intermediate-band filters 
(Fig. 1) covering a wavelength range from 3000 to 10 000 A, which are specifically designed for 
avoiding night sky emission lines. A detailed description of the BATC photometric system can 
be found in Fan et al. (1996) and Zhou et al. (2003). 

Fig. 1: Transmission of filters. 

From August 1995 to March 2002, we totally accumulated about 37 hours exposure for 
Abell 168 with 13 filters. More information about our observations is given in Table 1. The raw 
data are processed with an automatic data reduction software named PIPELINE I (Fan et al. 
1996) which includes the bias subtraction and the dome flat-field correction. The technique of 
integral pixel shifting is used in the image combination during which the cosmic rays and bad 
pixels are corrected by comparing multiple images. The mean error of position calibration is 
~ 0".5 by taking the Guide Star Catalog (CSC) as standard. By using the PIPELINE 11 software 
(Zhou et al. 2003), developed on the basis of the DAOPHOT kernel, photometries are performed 
for all the objects detected. As a result, the photometries with the Point Spread Function (PSF) 
model and the different apertures are obtained for 4798 objects detected in at least 3 bands. 
The flux calibration of the SEDs is performed by using the Oke-Gunn standard stars which 
were observed during photometric nights. The detailed information about calibration can be 
found in Zhou et al. (2001). Because we have no calibration image for g fllter, we instead 
perform a model calibration that is developed specially for the large-fleld photometric system 
by Zhou et al. (1999). 



Table. 1: The statistics of observations. 
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2.2. Combining the BATC and the SDSS SEDs of galaxies 

Abell 168 has been observed by SDSS with five broad bands, namely u' , g' , r', i' , z' (Fig. 1) , 
and the photometric data are distributed in the SDSS Early Data Release (EDR, see Stoughton 
et al. 2002). The method of combining the SDSS photometric data and the BATC SEDs is 
explored recently by Yuan et al. (2003), and this work shows that the combined SEDs could 
lead to a more accurate estimate of photometric redshifts. We cross identify the BATC source 
list with the SDSS photometric catalog in the same region, and a catalog containing 2580 
galaxies and 2218 stars is achieved. The object classifications are taken from the result of 
SDSS photometries. For the photometries of BATC, wc adopt the PSF magnitudes for stars, 
and take a fixed aperture with a radius of 4 pixels (i.e., rap ~ 6". 8) for galaxies. In order to 
combine the photometric results of the BATC and the SDSS, an aperture-correction should be 
applied to the SDSS model magnitudes (mmodei) by 

C'^'^ 2nrl{r)dr 

Am = map - mmodei = -2.5 log -joot; — w \j ' (1) 

Jq 27rr/(r)dr 

where map is the aperture magnitude, /(r) is the preferred profile function of surface intensity 
(e.g., the power law or de Vaucouleurs r^/^ model). The model magnitude in each band and 
its corresponding parameters that are used to quantify the preferred brightness profile can be 
found in the SDSS data. 

To observe the possible systematic offset between the BATC and the SDSS photometric 
system, we derive the zero point difference of SED (ZPDS) for each object by calculating the 
difference between the SDSS r' magnitude and the interpolated BATC magnitude at 6I8OA. 
Fig. 2 shows the distribution of ZPDSs for all the objects detected. For stars, the magnitude 
offset between these two photometric systems approximates to zero, indicating that the BATC 
and SDSS photometric systems are in good agreement. For galaxies, a more widely spread 
distribution of ZPDSs is found. This might be the result from the errors in photometries and 
the uncertainties in model fittings. For adapting to BATC multicolor photometric system, 
the SDSS SEDs are corrected by adding their ZPDS. Finally, we obtain the combined SEDs, 
including 13 BATC bands plus 5 SDSS bands, for 1553 galaxies brighter than r' — 20 mag. 

Fig. 2: Zero point distribution. 



2.3. Photometric Redshift and Membership Determination 

The photometric redshift technique was originally developed for detecting the high- 2; ob- 
jects based on the broad-band photometries (Pello et al. 1999; Bolzonella, Miralles, & Pello 
2000). The 15- intermediate-band SEDs obtained by the BATC photometric system can be re- 
garded as the rough spectrum. The advantages of the BATC data in measuring the photometric 
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redshift of some faint galaxies can be expected. With the BATC SEDs (Yuan et ah 2001, Xia 
et al. 2002) and the combined SEDs of the BATC and the SDSS (Yuan et al. 2003), the pho- 
tometric redshifts of faint galaxies can be estimated with a certain accuracy (A2;phot ~ 0.02), 
which allows us to isolate the member galaxies for some nearby galaxy clusters. 

The technique of photometric redshift has been applied to the 18-band SEDs with the help 
of the Hyperz code (Bolzonella, Miralles, & Pello 2000). The red-shifted spectral templates 
of normal galaxies are used to fit the observed SEDs, and the reddening law is set to be a 
free parameter during the fitting. The photometric redshift for each galaxy is searched in a 
redshift range from 0.0 to 1.0, with a searching step of 0.005. The histogram of Zp^^t for all 
galaxies is plotted in Fig. 3, and most of them are located at 2; < 0.5. 195 galaxies with known 
spectroscopic redshifts ^^gpec < 0.25 are selected to derive the uncertainty of photometric redshift 
which is found to be ~ 0.024 (see Fig. 4). 

Fig. 3: The histogram of photometric redshift. 
Fig. 4: A comparison of photometric redshifts and spectroscopic redshifts. 

According to the known spectroscopic redshifts, there are 121 galaxies distributed in a 
range from 0.04 to 0.05. We take them as the member galaxies of Abell 168. About 90% 
(109/121) galaxies have photometric redshifts within the range [0.02,0.07], thus we take this 
range as the criterion for the membership determination based on the photometric redshift. As a 
result, 255 galaxies are selected as the member candidates. Together with the spectroscopically 
confirmed members, we totally obtain an enlarged sample of 376 member galaxies. It should be 
noted that the uncertainty of photometric redshift is rather large compared with the intrinsic 
dispersion of Abell 168, but it still can be used to select the member galaxy candidates with 
high efficiency, especially for the faint galaxies. 

3. SUBSTRUCTURES 

3.1. Spacial Distribution of Members and X-ray Map 

Fig. 5 gives the spatial distribution of 376 member galaxies (small circles), with its corre- 
sponding contour map (thin lines) of the number density superposed. In the figure, D indicates 
the density peak located at 01^15°^06'^ +00°19'12" (J2000.0). The spatial distribution is elon- 
gated along the NW-SE direction with a position angle of ~ 30°. A possible group (A168WG), 
indicated by E in the figure, can be recognized to the west of the cluster core. The group is close 
to the high density region and seems to be falling down to the center, increasing the galaxy 
density there. A similar spatial distribution was earher derived by Kriessler & Beers (1997), 
but their sample of 106 galaxies did not show the core structures and the group clearly. We 
take the density peak as the center, and then derive the radial number density profile in Fig. 6. 
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It is clear that most of the galaxies concentrate within a range of about 12' (~ 0.7 Mpc). 
Fig. 5: The nicnil:)cr distribution and the X-ray map. 

Fig. 6: The number density profile. 

The X-ray image in 0.2 — 3.5 keV band, observed by IPC of Einstein Observatory, is 
overlayed on the spatial distribution in Fig. 5 (thick lines). It is interesting that the X-ray map 
and the member distribution coincide well with each other both in their spatial locations and 
their shapes. As mentioned above (§1), the X-ray emission of Abell 168 has two peaks, marked 
by A and B in Fig. 5. The higher X-ray peak, point A, is located at 01^14"^57" +00°25'04" 
(J2000.0) with a distance of ~ 6' (i.e., ~ 320 kpc) to the density peak. In the figure, C points 
to the cD galaxy UGC 00797. Point A is very close to the cD, which agrees with the traditional 
experience that D/cD galaxies are located around the peak of the X-ray emission (Jones et al. 
1979; Forman & Jones 1982). However, this cD galaxy significantly deviates from the local 
density peak of Abell 168, which seems inconsistent with the idea that a D/cD appears to 
be associated with the center of its host cluster (Beers & Geller 1983; Oegerle & Hill 2001). 
On the other hand, the lower X-ray peak, point B, is surrounded by the clustering bright 
galaxies. From the shape of the X-ray emission, the hot gas between the two peaks seems to 
be pushed away from the center by the head-on collision of the two clumps. Combining all 
above characteristics, it seems that there are some substructures in Abell 168. Ulmer, Wirth 
& Kowalski (1992) pointed out that Abell 168 was possibly formed by the collision of two 
subclusters. If this is true, it is reasonable to assume the collision occurs along the elongated 
direction. To explore the substructures in detail, we artificially divide the main region of 
Abell 168 (indicated by a 0.4° x 0.6° box in Fig. 5) into two equal parts. One is centered at the 
higher X-ray peak with 79 galaxies. We call it north subcluster (A168N). Another is located 
at the lower X-ray peak with 97 galaxies, named south subcluster (A168S). In the following 
sections we try to find other evidences to support that they are indeed physically different parts 
in Abell 168, in other words, they are really two subclusters, and we will show that they seem 
to be under a merging process. 



3.2. Dynamics of Substructures 

Provided Abell 168 is formed by merging of two subclusters, the dynamics remaining 
after the collision should be presented by testing the velocity distribution. On the basis of 
the spectroscopically confirmed member galaxies, we are allowed to study the dynamics of 
Abell 168 in detail. The statistical parameters of A168N, A168S, A168WG and A168 are listed 
in Table 2. The mean redshift (z), the mean velocity (vj.) and the velocity dispersion (dr) are 
all calculated by the biweight statistic in ROSTAT software (Beers, Flynn, & Gebhardt 1990), 
except the velocity dispersion of A168WG is taken from the gapping statistic. The asymmetry 
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index (AI) and the tail index (TI) of the velocity distributions are also provided by ROSTAT. 
Fig. 7 shows the velocity distribution of these two subclusters, respectively. The distribution 
of A168N is more concentrated than that of A168S. It seems that, for A168S, there is a gap at 
~ 13 000 km s~^, which can be clearly seen from a stripe density plot in Fig. 8. According to 
the value of TI (see Bird & Beers 1993 for detail), A168N possesses a hght-tailed distribution, 
while the distribution of A168S is near Gaussian. This suggests that these two distributions 
are intrinsically different. Moreover, the classical K-S test also shows the velocity distributions 
of these two subclusters are from the different parent distributions with 94% probability. 



Fig. 7: V, distribution of A168N and A168S. 



Fig. 8: Stripe density plot. 



Although it appears the obvious interaction between these two subcluster according to 
the spatial distribution, one should keep in mind this is a projection effect. So it is necessary 
to estimate the probability that they are gravitationally bound, in order to confirm they are 
really two interacting subclusters. We take a linear two-body model for the system as Beers, 
Geller, & Huchra (1982) did in analyzing the cluster Abell 98. Under this model, the relative 
velocity V and the separation R between two subclusters can be converted to the observational 
quantities: the relative radial velocity = l^sina and the projection separation Rp = Rcosa. 
The projection angle a (g [0,7r/2]) is defined as the the angle between the line connecting the 
two subclusters and the plane of sky. The energy relation for gravitational binding is taken 
from the Newtonian criterion 

V'^R 

< sin^acosa, (2) 

where G is the gravitational constant. M is the system mass, equal to the total mass of the two 
subclusters. The difference of mean velocity between A168N and A168S can be easily figured 
out Vj- = 264±142 km s^^. We take Rp = 0.584 Mpc, the distance between the two X-ray peaks 
which roughly represent the centers of these two subclusters. Following the method described in 
Beers, Geller, & Huchra (1982), the mass and the gravitational scale lengths of A168N, A168S 
and A168 are estimated and hsted in Table 2. The probability of gravitational binding can be 
obtained by the formula 

1 r+oo 

P-jI p{Vr)p{a\V,)dV,, (3) 

where A is re-normalization factor. p{Vr:) is the probability distribution of the velocity which 
is assumed to be a Gaussian. p{a\Vr), which can be obtained from Eq. 2, is the probability 
of the valid a for gravitational binding at a given V^-. Integration is over all the appropriate 
values of velocity. As the results, we find these two subclusters are gravitationally bound with 
a probability of 92%. 

Another efficient statistic method to test the substriictures of the radial velocity distri- 
bution was recently introduced by GoUess & Dunn (1996), who keep the spirit of Dressier- 
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Shectman test (Dressier & Shectman 1988) and employ the standard K-S test to compare the 
velocity distribution of the local group, n nearest neighbors, and the whole velocity distribution. 
The statistic is defined as 

N 

Kr^^Yl [^KS(^ > ^obs)], (4) 

i=l 

where is the number of all member galaxies of a cluster, and D is the statistic in standard 
K-S test. For each galaxy, we compute the probability, Pks{D > -Dobs), of its local group via 
Monte Carlo simulation by randomly shuffling velocities every time. This probability gives the 
an indication of how significant the observed £)obs is by comparing with the simulations. Kn is 
an indication that the local velocity is different from the overall distribution. 

We take the group size of n = 10. By involving 10^ simulations, we find that Kiq = 0.000%, 
almost no simulated D greater than the observed D'obs- A significant deviation of local velocity 
distribution from the overall distribution can be found in Fig. 9. The bubble radius for each 
galaxy is proportional to —log [Pks{D > -Dobs)]- Consistent with the result of spectroscopic 
redshifts, no further substructure was found in the central region, because the photometric 
redshift is not precise enough to distinguish the small difference, such as ~ 260 km s~^, in 
velocity. The result possibly suggests that the two subclusters are approaching each other 
along the plane perpendicular to our line of sight. The obvious bubble concentration at the 
west possibly supports the existence of the west group discussed in the section §3.1. 

Fig. 9: Bubble plot. 

4. LUMINOSITY FUNCTIONS 

4.1. Double Schechter Function 

The luminosity function is a key diagnostics for clusters of galaxies because it is tightly 
related to the dynamical evolution and the merging history of galaxy clusters. In the past 
decades, the LFs of GCs are well described by the famous Schechter function, which was 
originally proposed by Schechter (1976): 

(j){L)dL = (P* {L/L* f exp(-L/L*) d{L/L*), (5) 

where (f)*,L*,a are the normalization parameter, the characteristic luminosity and the faint 
slope parameter, respectively. 

Recently, the deep observation of the Coma cluster reveals their LFs possess an enhanced 
faint tail (e.g., Trentham & TuUy 2002). It is widely perceived that a single Schechter function 
gives a poor representation of the data when the LFs extend to the fainter end (Biviano et 
al. 1995; Durret, Gerbal, Lobo & Pichon 1999; Trentham & TuUy 2002). The segregation in 
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LFs of GCs suggests that there are at least two populations of the galaxies in clusters (Durret, 
Adami & Lobo 2002). For example, Ferguson & Sandage (1991) take a form of Gaussian part 
plus Schechter function in order to probe the two populations, the giants and the dwarfs. 

The Schechter function works well for the bright end as shown in the previous literatures, 
such as Dressier (1978), Paolillo et al. (2001). However, the present extending LFs require 
a variation to describe the enhanced tail rather than a single Schechter function. Hence we 
propose an acceptable double Schechter function (DSF hereafter). 

N{M) = <^*[lO-0.4(M-M*)]aoglO--(---o) 

+ 0*[lO-O-4(M-Mr)]a,glO-o-(---r)_ (g) 

Each term has the same form as the Schechter function in Eq. 5. M* {i = 0, 1) is the charac- 
teristic absolute magnitude with an equivalent characteristic luminosity; (pljCHi, {i — 0, 1) have 
the same meanings as in Eq. 5. 



4.2. Luminosity Functions 

We use DSF to fit the observed LFs by means of the minimization. The results are 
hsted in Table 3 and plotted in Fig. 10. A168S shows a decay tail while A168N shows an 
increasing one, suggesting that the galaxy contents of these two regions are quite different. We 

also compute the ratio of bright galaxies (M,./ < —18.5) to the faint [Mr' > —18.5) for these 
two subclusters. This ratio for A168N is found to be 32/47=0.68. However, for A168S, the 
ratio is 56/41=1.37. It is obvious that the bright galaxies are more abundant in A168S while 
the faint galaxies prefer to occur in A168N. This can be interpreted by the cannibalism model 
(Hausman & Ostriker 1978). That is to say: the cDs are likely to be formed by accreting the 
surrounding massive galaxies, which naturally results in the deficit of bright galaxies in the 
most dynamically evolved clusters (Lugger 1986). Hence, the formation and evolution of these 
two subclusters might be different. A168N, associated with the cD galaxy UGC 00797, is hkely 
to be more evolved than A168S. To confirm our dividing of A168N and A168S is not accidental. 
We try other methods of dividing the core region. For example, we divide the core region into 
an east and a west part equally along the elongated direction of the spatial distribution. In this 
case, the clear contrast in the LFs disappears, which supports the dividing of A168N and A168S 
is a special elioiee. arid tliey are tlie Uxo sul)elusters under tlie merging e\"orit of Al)ell 168. 

Table. 2: The fitting results of LFs of Abell 168. 
Fig. lOrThe observed luminosity functions are fitted by DSF. 

Panel (c) in Fig. 10 shows the LF of the central region (A168N+S) of the cluster. As shown 
in previous works of Dressier (1978) and Oegerle, Ernst, & Hoessel (1986), the LF has a quite 
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flat tail which is, however, an average of LFs of A168N and A168S according to the present 
work. As for the whole field, the observed LF and its fitting model are shown in panel (d) of 
Fig. 10. For the overall LF, DSF gives a good description to the obvious enhanced faint tail at 
Mj-i > —18.5 mag. Compared with the LF of the central region, the overall LF suggests that 
the faint galaxies (possible dwarfs) tend to be located in the outer region of the cluster. This 
perhaps support the population-density relation: dwarfs are more common in lower density 
environments (Phillipps, Driver, Couch, & Smith 1998). 

5. SUMMARY 

We accumulated about 37 hours of high quality observations of Abell 168 with 13 interme- 
diate filters of BATC photometric system. In the EDR of SDSS, we find that the photometry 
and the spectroscopic redshift measurement are completed in the same field. After a cross- 
identification with SDSS, we combine the photometric data of the objects detected in two 
systems getting the 18-band SEDs, which include 13 BATC bands and 5 SDSS bands. A 
complete sample of 1553 galaxies brighter than r' — 20 mag is achieved for the subsequent 
study. 

By using the technique of photometric redshift, we get an enlarged sample of 376 galaxies, 
including 121 spectroscopically identified member galaxies and 255 newly selected member 
candidates. Spatially, the member galaxies show an elongated distribution along the NW-SE 
direction with a position angle of ~ 30°, which is consistent with the spatial characteristic of 
the X-ray image obtained by Einstein Observatory. Moreover, the two X-ray peaks and the 
deviation of the cD galaxy UGC 00797 from the local density peak support the existence of 
substructures in Abell 168. According to the spatial characteristics, we artificially divide the 
core region into two parts, A168N and A168S. By employing ROSTAT software, a detailed 
investigation of the dynamics of these two subclusters is performed. A relative radial velocity 
of 264 ± 142 km s~^was found between them. Their velocity distributions are intrinsically 
different and should be from the different parent distributions. Via building a linear two-body 
model, they are found to be gravitationally bound with a probability of 92%. Investigation of 
LFs by the DSF shows an obvious difference in the galaxy content of these two subclusters. 
The A168S is dominated by bright galaxies while A168N has larger number of faint galaxies. 
A168N is a sub-system associated with the cD galaxy and likely to be a more evolved system 
than A168S. All these indications are consistent and support that Abell 168 is hkely to be a 
head-on merging system at its early stage of the process, and that the collision plausibly occurs 
along the direction perpendicular to our line of sight. 
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Table 1: Parameters of the BATC filters and the statistics of observations of Abell 168. 



No. 


Filter 


^eff 


FWHM 


Exposure 




Seeing" 








(A) 


(A) 


(hours) 




(arcsec) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


c 


4210 


309 


3.67 


11 


5.08 


1 


2 


d 


4546 


332 


5.33 


16 


4.09 


1 


3 


e 


4872 


374 


4.00 


12 


3.94 


2 


4 


f 


5250 


344 


3.47 


11 


4.18 


4 


5 


9 


5785 


289 


1.33 


4 


3.81 





6 


h 


6075 


308 


1.00 


3 


4.80 


3 


7 


i 


6710 


491 


1.08 


4 


3.80 


3 


8 


3 


7010 


238 


2.00 


6 


3.94 


1 


9 


k 


7530 


192 


2.00 


6 


4.49 


3 


10 


m 


8000 


255 


5.67 


17 


3.94 


1 


11 


n 


8510 


167 


2.67 


8 


4.81 


2 


12 





9170 


247 


1.33 


4 


6.29 


1 


13 


P 


9720 


275 


3.40 


11 


5.03 


2 



"^Number of images of each color. 
''The seeing is of the combined images. 
■^Number of the calibration images. 
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Table 2: Parameters of Abell 168.'^ 



Parameter 


A168N 


A168S 


A168WG 


A168 




79 


97 


19 


376 


-^gal, spec 


26 


52 


6 


121 


Z 


0.0446 ± 0.0004 


0.0455 ± 0.0003 


0.0466 ± 0.0001 


0.0451 ± 0.0002 


(km s"""^) 


13065 ± 113 


13329 ± 86 


13633 ± 35 


13206 ± 50 


C7r (km s^^) 


564 ± 90 


613 ± 56 


491 ± 189 


554 ± 34 


(l/rp)-i (Mpc)^ 


0.30 


0.30 




0.77 


M (IQi^M©) 


2.1 ± 0.7 


2.5 ± 0.5 




5.2 ± 0.6 


AI^ 


-0.263 


-0.392 


-0.327 


-0.093 


TP 


1.458 


0.958 


0.000 


0.901 



'I'The errors arc computed at 68% confidence level. 

'^The number of galaxies including the spectroscopically confirmed members and the photometric redshift selected 
candidates. 

^The number of galaxies only including the spectroscopically confirmed members. 
■^The gravitational scale length. 
''Asymmetry index. 
®Tail index. 
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Table 3: The fitting results of LFs of Abell 168. ^ 









Mo* 






Ml* 




• 

•A- mm 


A168N 


14.9 


+3.4 
-2.9 


-21 


0.041H? 


07 r+20.2 


-16.911°:^^ 


40+°-^* 
'-'•^'-'-1.00 


1.6 


A168S 


54.4 


+8.5 
-7.4 


-20.47t°iJ 


0.36t°;0« 


0.00 






3.4 


A168N+S 


50.3 


+5.6 
-5.2 




•^•^0-0.04 


0.00 






4.1 


Total 


95.8 


+10. 
-9.4 


-20.77l°:23 


n 1 7+0.05 
^•■'■'-0.05 


138.9t24-5 


-17.48l°i^ 


^•J^^-0.26 


3.6 



^'The errors are computed at 68% confidence level. 
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4000 6000 8000 10^ 
Wavelength (A) 



Fig. 1. — The transmission curves of filters in the BATC (sohd fines) and the SDSS (dashed 
lines) photometric system. The names of the filters are marked at the top of each curve. 
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Fig. 2. — The histogram of the zero point differences of individual SEDs in the BATC and the 
SDSS system for all galaxies and stars. 



-19- 



o 

CO 



o 
u 



o 

C\2 




0.2 0.4 0.6 0.8 0.02 0.04 0.06 0.08 

Photometric redshift 

Fig. 3. — Panel (a) shows the histogram of photometric redshift for 1553 galaxies. The distri- 
bution from 0.0 to 0.1 is shown enlarged in panel (b). In both figures, the solid lines represent 
photometric redshifts of all galaxies, the dashed lines indicate photometric rcdshifts of the 
galaxies without spectroscopic redshift. The range of [0.02, 0.07] marked by two arrows is 
taken as the criterion for member candidates. 
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Spectroscopic redshift 



Fig. 4. — A comparison of photometric redshifts and spectroscopic redshifts for 195 known 
galaxies. The dashed hnes mark out the two-sigma region of the photometric redshifts. 



01^17"^00^ qqs i6°^00^ 30^ 15™00^ 30^ 14°^00^ 30^ 



RA (J2000.0) 



Fig. 5. — Member distribution (small circles) of Abell 168, with corresponding contour map 
superimposed (thin lines). Point A and point B are the two X-ray peaks. C points to the cD 
galaxy UGC 00797. D points to the number density peak. E points to the possible west group 
indicated by a dashed hne ellipse. The box (0.4° x 0.6°) indicates the central region of the 
cluster. The dotted-dashed line divides the north subcluster (A168N) and the south subcluster 
(A168S) equally (see §3 for detailed analysis). 
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Fig. 6. — The radial number density profile. The bottom abscissa shows the distance in unit 
of Mpc. The top abscissa is in unit of arcmin. 
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Fig. 7. — The velocity distributions of the spectroscopically confirmed member galaxies 
A168N and A168S are presented by the solid line and the dashed line, respectively. 
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Fig. 8. — For the spectroscopically confirmed members in these two subclusters, the stripe 
density plot of velocity is plotted above the spatial distribution. 
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RA (J2000.0) 

Fig. 9. — The bubble plot shows the deviation of the local velocity distribution for groups of 
10 nearest neighbors from the overall velocity distribution. By involving 10^ simulations the 
significance of existence of substructures reaches P{kiq > k'^^) = 0.000%. 



-26- 



CV2 - 



O 



CV2 - 



T — I — ' — I — ' — r 
(a) A168N 



:(c) A168N+S 




5*5 =- 



T — I — ' — I — ' — I — ' — r 
(b) A168S 




(d) Total 




-22 -20 -18 -16 -22 -20 -18 -16 



Absolute magnitude {M ) 



Fig. 10. — The observed luminosity functions are fitted by DSF. The dashed lines and the 
dotted lines indicate the components of DSF and the sohd fines are the sum of them. Panel (a) 
and (b) are of the north subcluster and south subcluster, respectively. The LF of the central 
region is plotted in (c). In panel (d) we show the LF of all 376 members in the whole field. For 
A168S and A168N+S (equal to the central region of Abell 168), the single Schechter function 
gives a better representation than DSF. The detailed discussion on LFs is made in §4.2. 



